Polar Kerr Effect and Time Reversal Symmetry Breaking in Bilayer Graphene 
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The unique sensitivity of optical response to different types of symmetry breaking can be used 
to detect and identify spontaneously ordered many-body states in bilayer graphene. We predict 
a strong response at optical frequencies, sensitive to electronic phenomena at low energies, which 
arises because of nonzero inter-band matrix elements of the electric current operator. In particular, 
the polar Kerr rotation and reflection anisotropy provide fingerprints of the quantum anomalous 
Hall state and the nematic state, characterized by spontaneously broken time reversal symmetry 
and lattice rotation symmetry, respectively. These optical signatures, which undergo a resonant 
enhancement in the near-infrared regime, lie well within reach of existing experimental techniques. 



Optical experiments have been successfully used to 
probe diverse electronic phenomena in graphene For 
bilayer graphene (BLGh physical properties such as the 
gate tunable bandgap [2|, |3j , the band structure param- 
eters and the electron phonon coupling 0, E| were 
investigated with the help of infrared and optical spec- 
troscopy. These techniques have also been used to probe 
interaction effects such as band renormalization [3, Hoj 



and exciton formation 11, 121. However, there has not 



yet been any effort to apply optical methods to the in- 
vestigation of strongly correlated states, which are ex- 
pected to form in BLG at low energies [131421 1 . This can 
be partly due to the low characteristic energy scales for 
these symmetry breaking states, estimated to be of order 
1 me V 14( , which lie far outside the range of character- 
istic energies probed in optical experiments. 

In this Letter, we point out that the problem of en- 
ergy scales is offset by the unique sensitivity of optical 
response to broken symmetries. This, along with sev- 
eral other features, makes these methods ideally suited to 
the investigation of the interacting ground state of BLG. 
The possible broken symmetries are expected to mani- 
fest themselves through characteristic transport proper- 
ties such as a non-zero Hall response or anisotropy in lon- 
gitudinal conductance 13h21|. Detecting these effects in 
transport experiments requires fabrication of samples of 
BLG with at least four contacts, which proves challenging 
in suspended BLG currently used in these experiments. 
However, optical experiments allow us to measure the 
AC conductivity in a contact-free manner. As we discuss 
below, the AC conductivity shows distinctive signatures 
of broken symmetry just like the DC conductivity. These 
signatures are strong due to nonzero inter-band matrix 
element of the current operator for BLG. Thus the opti- 
cal response, which features additional resonant enhance- 
ment in the near-infrared regime, can be used to directly 
probe spontaneously broken symmetries in BLG. 

A large number of possible interacting phases have 
been proposed for BLG 13l42l|. Recent compressibil- 
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FIG. 1: Kerr angle (in units of fine structure constant a = 
e 2 /he) as a function of photon energy for BLG in the QAH 
phase. Note the resonant enhancement near Eo = 0.4 eV, 
arising from direct transitions to the higher BLG bands, 
Eq. (|12[) . Inset: Schematic band structure of BLG near the 
K point, for the QAH phase. The Kerr response arises from 
transitions 1' — > 2 and 1 — »• 2', involving states in the bands 
1 and 2 which are affected by broken TRS. 



ity and transport experiments on charge neutral, sus- 
pended, double gated bilayer graphene [22t - 424j j appear to 
confirm the prediction of a non-trivial interacting ground 



state. The experimental data was argued [24{ to be con- 
sistent with only two of the proposed phases: the Quan- 
tum Anomalous Hall phase (QAH) predicted in (la. Il7|. 
and the nematic phase predicted in flil - Eoj ] . Both these 
phases are uniquely interesting phases. The QAH phase 
spontaneously breaks time reversal symmetry (TRS) and 
exhibits quantum Hall effect at zero magnetic field, while 
the nematic state involves a distortion of the Dirac band- 
structure that spontaneously breaks the exact rotational 
symmetry of the lattice. If either of these phases is con- 
firmed in BLG, it would fulfill a long quest for an experi- 
mental realization of a QAH instability [2|| (QAH phase) 
or a Pomeranchuk instability [26], [27[ (nematic phase). 

Optical methods are ideally suited to identifying the 
ground state of BLG. The polar Kerr effect, wherein 
linearly polarized light has its polarization axis rotated 
upon reflection, is a well known optical probe of the Hall 
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conductivity. It has been used to probe quantum Hall 
states [28j, and more recently has been applied to topo- 
logical insulator thin hims in the vicinity of a ferromagnet 
23 . 30l j . and to p + ip superconductors 31 1. 

As we shall see, the QAH phase exhibits an AC Hall 
conductance in addition to the quantized DC Hall con- 
ductance, and thus the Kerr effect offers a direct test 
of the QAH scenario for BLG. Our analysis of optical 
response, taking into account transitions between four 
BLG bands, reveals a resonant enhancement of the AC 
Hall conductivity (see FigflJ. This resonant enhance- 
ment occurs because the microscopic current operator 
has inter-band matrix elements (FigfT] inset) correspond- 
ing to transitions from the low energy bands 1 , 2 to the 
high energy bands V , 2'. The resulting, resonantly en- 
hanced Kerr rotation is plotted in FigJTJ The predicted 
effect is many orders of magnitude larger than that ob- 
served in p-wave superconducting materials |3l( , and lies 
well within reach of existing experimental techniques. 

Optical methods can be used to probe domain forma- 
tion expected to occur in the TRS breaking QAH phase. 
Since different domains will produce a Kerr rotation of 
opposite sign, the spatial domain structure can be di- 
rectly imaged in optical experiments - a significant ad- 
vantage over transport experiments, which can only mea- 
sure the net effect of all domains. For a non-focused op- 
tical experiment, the effect of random domains will be to 
reduce the total Kerr angle by a factor \/Wo, where No 
is the number of domains. 

While the Kerr rotation allows to test for TRS break- 
ing, anisotropy in reflection allows to test for rotation 
symmetry breaking. As we discuss below, this leads to a 
characteristic dependence of the reflection amplitude on 
the polarization angle of incident light which offers a way 
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to test the nematic scenario for BLG 

Finally, we note that spontaneous symmetry breaking 
is only expected to occur below a critical temperature, 
estimated to be of order 1 - 10 K [13, [Hj]. The optical 
signatures of interacting states will thus show a strong 
temperature dependence, and will vanish entirely above 
a critical temperature. This provides a way to distinguish 
spontaneously broken symmetries from explicitly symme- 
try breaking effects (e.g. magnetic impurities), which will 
not show any comparable temperature dependence. 

Electron properties of a clean BLG are governed by a 
four-band Hamiltonian written for the four component 
wavefunction ip — (ipi, V>2, ip3, ipi), describing electron 
wavefunction on the sublattices A, B and A', B' on the 
two layers: 
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-ipei + e-ipej^ where f Q ~ 3.1 eV is 



the upper and lower bands. The quantity t p vanishes at 
the K and K' points, behaving as vp+ near point K and 
as — vp— near point K', where p± = p x ± ip y . 

The Hamiltonian (fT]) features four bands with energies 
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Near the points K and K', this gives two massless Dirac 
bands e\^{v) that cross quadratically at zero energy, and 
two high-energy bands £3,4 (p) ~ ±Eq. The dispersion 
near K and K' can be obtained by expanding in small 
tp/E Q , giving £i j2 



= ±|ip|7#o = ±v 2 p 2 /E , e v ,2> 



±{E + v*p 2 /E ). 

We now consider the effect of interactions. Interactions 



can open a bulk bandgap between bands 1 and 2 13Ml7| , 
resulting in a bandstructure of the form FigfTJinset). One 
particularly interesting gapped state is the QAH state, 
UMUM, the mean field Hamiltonian of which we present 
below. To exhibit more clearly the block structure we 
reorder basis vectors by interchanging the components 
ip2 and "04 ■ In this representation, we obtain 
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where A is the order parameter describing gap open- 
ing at the K and K' points. The other possible gapped 
states 13 -lH have a similar mean field Hamiltonian, 
but the sign of A is distributed differently among the 
spins and valleys. We note that under time reversal, 
H K {A)@H K ,(-A) -> H K ,(A)®H K (-A), so this phase 
breaks TRS. In consequence, the QAH state can exhibit 
a non- vanishing Hall conductance at zero magnetic field. 
However, the gap preserves the isotropy of the bandstruc- 
ture. Thus, the QAH state must exhibit isotropic longi- 
tudinal conductivity. 

Next, we discuss the relation between the Hall response 
in the QAH phase and the Kerr rotation. We consider 
an experimental setup where light is incident normally 
on a BLG sheet that is placed on a substrate with re- 
fractive index n, which is taken to be real (complex case 



considered in [33j). If the BLG sheet has a non- vanishing 



the hopping amplitude, and Eq is bandgap parameter for 



Hall conductance, then incident linearly polarized light 
will be reflected as elliptically polarized light, with the 
major axis of the ellipse rotated with respect to the in- 
cident polarization by the Kerr angle 9k- The standard 
formula relating the Kerr angle to the Hall conductance 
is 9k ~ Im ip'yx) [13] ■ However, this formula is derived 
for light incident on a conducting half space, whereas we 
are considering a BLG sheet that is much thinner than 
the optical wavelength. For this case, the relationship 
between Hall conductivity and Kerr angle must be cal- 
culated afresh, by solving the Maxwell equations on two 
sides of the BLG sheet and matching solutions at the 
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boundary. We obtain [33j 
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We now calculate the magnitude of the Kerr rotation, by 
evaluating the conductivity. The AC conductivity can be 
written using the Kubo formula as 
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where i and j are band indices, = n(si^ p ) is a Fermi 
function, and the sum over momenta p stands for an 
integral. The velocity operators V a are defined as V a — 
dH(p)/dp a , and 7 describes the excited state lifetime. 

We focus on the contributions which correspond to op- 
tical interband transitions between the massless low en- 
ergy bands (i=l,2), and the high energy bands (i=l',2'), 
which are separated from the low energy bands by the en- 
ergy E . We focus on these transitions because they are 
of resonant character at a frequency close to the band 
separation energy Eg, and hence dominate the optical 
response. We now note that {i,p\V a \j,p){j,p\Vp\i,p) = 
Ti(y a Ili p Vj3llj. p ), where ±T jP projects onto the state in 
band i with momentum p. Assuming we are at a tem- 
perature T < E /k B = 4000X, we obtain 
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where we used the relation £2' — £1 = £2 — £i' that follows 
from particle/hole symmetry of the Hamiltonian ([3]). 

We evaluate the expression ([6]) for p near point K with 
the help of the projectors 
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Here 2' project on the Bl and A2 sublattices (lower 
right corner of the Hamiltonian in ©), and a x acts on 
this subspace. Meanwhile, III 2 project on the Al and 
B2 sublattices (upper left corner of Eq.©), and h(p) 
is the effective two band Hamiltonian for the massless 
Dirac states, which has eigenvalues E(p) — ±\\h(p)\\. 
The trace over projectors takes the form 
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Here V Q denotes d/dp a . We now compute h(p) by using 
second order perturbation theory in vp/En, and obtain 



hx{p) = 



A 

v 2 p 2 _/E a 



v 2 p 2 + /E 
-A 



h K ,(A) = h* K (-A). 

(9) 

This result agrees with |34| . We substitute this two band 
Hamiltonian into Eq.© and obtain 



iG=gS - K 1 ~ iw) Vwt>)v - tCp) (10) 
+ I ( 1+ p5i) v ^ (p)v/(p) ' 

where we suppressed the terms arising from off-diagonal 
parts of h(p) — these terms give zero upon integration 
over d 2 p. Hence, we find g 12 y = -g 2 xy x = \m 2 ^^ y We 



substitute these results into Eqs.©,©, to obtain 
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where ./V = 4 is the number of spin/valley flavors, and 

^P = S 2 (p) ~ £i'(p). 

We now specialize to optical frequencies ui ^> A, 
and also assume 7 > A. We approximate by taking 
\\h{p)\\ « v 2 p 2 /E and O p « £ + 2v 2 p 2 /E a , and per- 
form the momentum integral in polar co-ordinates. The 
log divergence near p 2 = is cut by |A|, but there is 
no need for any high energy cutoff. In this manner, we 
obtain the Hall conductivity 
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There is also a contribution from 1 — »• 2 



which may be evaluated in the two band model [21 



transitions, 
II. 



This contribution, extrapolated to optical frequencies 
lj — E > A is of order (A| A\/uj 2 )e 2 /h. This is smaller 
than the contribution (| 12[) by a large factor 

fl4>l. (13, 

Thus, the Hall conductivity at optical frequencies is dom- 
inated by transitions to the higher bands, necessitating 
our four band analysis. 

From the result Eq. (fT2")) and the expression Eq.© we 
can extract the Kerr angle Ok- We take n — 1.5, which 
describes SiC-2 substrate, and take A s» 1CP 3 eV [1J|, and 
7 = 0.05 cV. The resulting Kerr angle as a function of 
frequency is plotted in FigJT] In optical experiments on 
cuprate materials, Kerr angles as small at 10 ~ 9 have been 
measured (3l| . The six orders of magnitude larger Kerr 
rotation in the QAH phase should thus be comfortably 
within reach of experiments. 
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The nematic state Tsj~20] is another interestini 



I or 



dered state proposed to explain the experiments [23j, [24 \ . 
This state is time-reversal invariant, featuring no Kerr ro- 
tation. Instead, it breaks rotation symmetry of graphene 
crystal lattice. The Hamiltonian for this state is 
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After reduction to the two low-energy bands, it becomes 
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(15) 

This Hamiltonian describes splitting of the quadratic 
band crossing into two linear band crossings. The ar- 
gument of the nematic order parameter A specifies the 
orientation of the nematic axis, which is defined as the 
line joining the two linear band crossings. The nematic 
axis makes an angle ip = —tt/2 + arg(A)/2 with respect 
to the p x axis. The nematic state manifestly breaks the 
approximate rotation invariance of the low energy band- 
structure, which manifests itself in an anisotropic longi- 
tudinal conductivity. Writing cr(6) = <tq + 8a{9), where 
is the angle with respect to the x axis, we obtain an 
expression for the reflection amplitude r{9), 



"(0) 



1 — n 
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For high frequencies lu 3> A, we calculate using the for- 
malism introduced above that 



h u> A 



cos(2(0-p)) (17) 



Again, this exceeds the anisotropy calculated in the two 
band model [21| by the large factor Eq. (fl"3|) . We note 
that trigonal warping of the BLG bandstructure arising 
from higher neighbor hopping can also lead to a reflec- 
tion anisotropy. However, these effects respect the three- 
fold rotation symmetry of the lattice. In contrast, the 
anisotropy resulting from formation of a nematic state 
exhibits a twofold rotation symmetry. The breaking of 
the exact lattice rotation symmetry can serve as diagnos- 
tic of the nematic state. 

To conclude, optical experiments can be used to probe 
broken symmetries in BLG by measuring the conductiv- 
ity in a contact free manner. The polar Kerr effect, by 
providing a means for measuring Hall conductivity, can 
be used to detect the QAH phase. TRS breaking gapped 
states that do not display a Hall conductance [17| can 
also be probed using the Kerr effect, although for these 
states the Kerr angle will be smaller than that for the 
QAH state by the small parameter d/X, where d = 3A 



is the BLG interlayer spacing and A is the wavelength 
of the light used in the experiment {36} . Nevertheless, 
this much weaker Kerr rotation will still be much larger 
than that measured in 31) . and will be within reach of 
experiments. Meanwhile, the nematic scenario for BLG 
may be probed by looking for an angle dependence of 
the reflection amplitude, which provides a direct test of 
broken rotational symmetry. 
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For an incident wave polarised along the x axis, E- m = 
x, we have E> = t xx x + t yx y, E< = (1 + r xx )x + r yx y, 
H> = n{t xx f-t yx ±), i?< = (1 - r xx )y + r yx x. Applying 
Ampere's law to the magnetic field leads to the continuity 
relations 



^ra 1 "F T x 



47T 



" : r ' .r.r — ^ {f^xx^xx H~ Oxytyx) 



^ityx ^yx 



Air 



{cfy X t X x "F &yytyx) (1^) 



Eliminating t using the continuity relations for electric 
field, E < — E>, we obtain the single matrix equation 



APPENDIX 

We consider light incident normally on a BLG sheet 
placed on a substrate with refractive index n = n' + in" . 
Incident and transmitted waves propagate in the +z di- 
rection, while the reflected wave propagates in the — z di- 
rection. The BLG sheet is taken to be in the z — plane, 
whereas the substrate occupies the halfspace z > 0. 
We calculate the reflection amplitudes for incident light 
that is linearly polarised along the x axis. The re- 
flected wave, E r = r xx x + r yx y, is linearly polarised 
if Im (r yx /r xx ) — 0, and elliptically polarised otherwise. 
The major axis of the polarisation is rotated with respect 
to the x axis by the Kerr angle 9k — Re (r yx /r xx ). 

We start with rewriting Maxwell's equations k x H = 
-^D,kxE=^Has 



kH x = n E y 



kH y — —n E x 



(18) 



kE x = —H,„ kE,. 
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Similar relations hold in the vacuum region with n re- 
placed by 1. At the interface z = we must match EM 
field amplitudes on both sides using continuity of the E 
field E < = Ey and the Ampere's law for the H field: 

Air Air 

(H>-H < ) x = -(*E) y (H > -H < ) y = --(aE) x . 
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This equation can be solved to obtain 
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We have denoted D = (n + 1 + ^o xx ) 2 - \ ^-a xy \ 2 w 
(n + l) 2 , and have assumed isotropy, so that a xx — <j yy . 
The Kerr angle is given by Re (r yx /r xx ), and thus takes 
the form 
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c[(n' 2 - n" 2 - l) 2 + 4n' 2 n" 2 ] 

where in the last line we have taken n = n' + in" with n' 

= 0, we 



and n" real. Now if we assume 



and 



obtain the formula quoted in the main text 
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